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Abstract 


A reverse slow shock and a forward slow shock were found in plasma 
and magnetic field data from Pioneer 6. The shocks were oblique an'> vc^ak 
(Mach «1.2) . 

Numerous non- shock- like discontinuities were found with B^f^B^ard 
^1^2* Assuming pressure balance across them, we obtained an average 
electron temperature * (1.8+,2)xlO^®K. It was observed that 
discontinuities with B^f^B^, ever occur. 
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I. Introduction 

Recently, a reverse perpendicular fast shock (Burlaga, 1970) and 
two forward slow shocks (Chao and Olbert, 1970) were identified in the 
solar wind. There is also some evidence for a reverse oblique fast 
shock (Chao and Blnsak, 1971). Reverse slow shocks have not hitherto 
been reported. 

We have scanned Pioneer 6 data in seach of such shocks 
with the aim of confirming and extending the above results. The 
search is described in Section II and the results are discussed Ir. 
Section III, Numerous discontinuities were found which are not shocks. 

From these we obtained an estimate of the average solar wind electron 
temperature which is independent of spacecraft sheath effects (SectionlV), 

The data used in this study are from Pioneer 6 which was launched 
on December 15, 1965 into a heliocentric orbit with apogee .98 AU and 
perigee .82 AU. A plot of its orbit during the early days of the flight, 
which turn out to be the times of interest, may be found in Ness (1966). 
Tlie magnetic field data were obtained by Ness and the plasma data by 
Bridge, Lazarus, and Davis. A discuss-lor. of their instruments may be 
found in Ness ^ (1966) and Lazarus ^ (1966), respectively. 
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II. Selection of Discontinuities 

The search for shocks was begun by looking for discontinuities 
in the magnetic field intensity B, since a jump in B is a necessary 
condition for all shocks except a parallel shock and since B is measured 
more frequently than the plasma parameters. A quantitative criterion 
for selecting changes in B was not used, since the variations in B are 
very complex. Instead, they were selected subjectively with an attempt 
to choose mostly step-like changes, the change being +207o in <1 min and 
the field intensity be ng nearly uniform over the. 5 min intervals before and 
after the change. Plasma data were not available for all of these 
discontinuities. The results that follow concern only those discontinuities 
for which both plasma and magnetic field data are available. 

To isolate the shocks, we first selected those discontinuities 
which had the correct signature (Colburn and Sonett, 1966), then 
eliminated those which clearly did not satisfy the Rankine-Hugoniot 
equations using the fitting scheme of Chao (1970), and finally eliminated 
those discontinuities which could conceivably satsify the pressure balance 
condition for tangential and rotational discontinuities. Discontinuities 
which were not shocks were considered spearately and were used to obtain 
the electron temperature as discussed in Section IV. 

The distribution of B for the discontinuities which were selected 
from Pioneer 6 data for this study is shown at the top of Figure 1. The 
absence of discontinuities near B * 0 is due to our selection criterion. 

The corresponding distributions of density changes and thermal speed 
changes are also shown in Figure 1. The changes- in n. were <107. for less 
than 5^ of the discontinuities'; This has implications concerning 
some results presented recently by Ivanov (1970). 
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Ivanov pointed out that In an anisotropic plasma, the hydromagnet Ic 
equations allow the exlstance of a kind of discontinuity with the 
following characteristics: 1) there Is a mass flux through the surface, 

2) the magnetic field Intensity changes across the surface while the 
density Is the same on both sides of the surface, 3) there Is a component 
of B normal to the surface, and 4) the components of B parallel to the 
surface are themselves parallel. He calls this a "rotational" discontinuity 
but this term Is misleading since the tangential component of ^ does not 
rotate. In fact, it is aao<re like a "shock" for vhicl aense that 

tt satisfies the coplanarity theorem and other shock equatimiis. 

Figure 1, shews that <57o of the discontinu ties which are identified 
in Pioneer 6 data satisfy the condition B^^ ^ B^, n^ * n^ necessary for 
Ivanov's discontinuities. Thus, Ivanov's discontinuities occur infrequently 


if at all. 
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IIL Slow Shocks 

For A flow shock, the following conditions must be satisfied. 

1) The measured 14 parameters, B^, Vg, n^ and n^ must 

satisfy the Rankine-Hugoniot equations. 

★ 

2) The computed Alfven Mach numbers, ahead and behind 

★ 

the shock must be both less than one, where is the component 
of^^'Che^ bulkt velocity normaiito the shoakofkoAtoandr^meaausad 
in the shock frame of reference, and is the Alfven speed 
based on the magnetic field component normal to the shock front. 

3) The Mach nuzibers based on the slow mode of magnetosonic waves, 

★ 

M , -V /V - , must be larger than 1 ahead of the shock and 

slow n slow 

less than 1 behind the shock; here V , is the speed of the 

slow 

slow mode of the magnetosonic wave in the direction of the 
shock normal. To compute one must know the total 

pressure of the plasma. Although this pressure is not measuie d 
directly, it can be es4;;.ipiated from the conservation equations 
for the normal momentum flux and the energy flux. 

4) Tlie computed electron temperatures ahead and behind the snock 
should be reasonable, i.e., on the order of 1.5x10^^. 

(Montgomery et al., 1970; Hundhausen and Montgomery, 1971). 
Besides showing that these conditions are satisfied, it is desirable 
to show that the events do not satisfy the pressure balance condition for 
tangential or rotational discontinuities. 

A Reverse Slow Shock. Figure 2 shows a discontinuity which has 
the signature of a reverse slow shock, namely, a discontinuous decrease 
in n and T and a simultaneous increase in V and B. It was seen at 
Pioneer 6 at 0109 UT on January 19, 1966. Using the procedure described 
by Chao (1970), a set of 14 parameters was obtained (^j^, Bg, Vg, n^^. 
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and see Table I) which 

•aclsfy the Rank Ine-Hugon lot equations for an Isotropic plasma and 

%fhlch are close to the observed average values for the discontinuity. 

We use here the RTN coordinates In which the R axis Is out from the 

sun and parallel to the sun-earth line, the T axis Is In the direction 

of the motion of the earth, the R~T plane is parallel to the ecliptic, and the 

N*»axls li' not^hvard and^f erpandlicular to the ecllptic^^ tliftke "th^9r%tical" 

vfcl«eet are* ‘fe'l veil* ifii> Table'^ 1 and arte ahowhi'byl the" sb^lid^ lines ite'^iigCiTe 2, They 

U9e <tot<>e, aoni^ared with the observed valnee^ showi^ln Figure' 2", The close 

agreement ahows that the observed parameters and changes are consistent 

with the Ranklne-Hugonlot equations. Using the RaDldnia»HugiDcilot v"! '* 

equations, it Is possible to solve for 3 of the 14 measured parameters 

as a function of the remaining 11 (Chao and Olbert, 1970). In particular, 

one can solve for as a function of (B, , B^, V,, n, , n_). Then, the 

measured parameters , n^ and are allowed to vary 

Independently within their uncertainties for computing the values. 

When the computed value of Is within the uncertainties of Its measured 

average value, then this computed value can be regarded as a prediction. 

Table 2 compares the ’’predicted" and observed values of The 

agrrement is good, c(See Figures 2, and cB 

The velocity of the shock surface relative to the spacecraft was 

found to be V *(212, 1 , 135)km/sec, corresponding to a speed 251/km/sec 
~ s 

A"' 

and a shock normal n*(0.85, 0.01, 0.54)^ rtiu i,«- lo i t* spactcraif 

The angle, 0_ a, between B. and n was found to be 42° as given in Table 
D • n *^1 

3. Relative to the shock, the normal component of the flow velocity 
decreased from 29 km/sec to 23 km/sec (See Table 3). The Mach numbers 
and ahead and behind the shock are given in Table 4. They 

do satisfy the conditions for a slow shock, l.e. the flow Is sub-alfvenlr 
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on both sides of the shock, It is super-slow ahead of the shock and it 
is sub'slov behind the shock. 

This discontinuity cannot be a tangential discontinuity with 

T «T , for the pressure balance condition gives a negative electron 
®1 ®2 

temperature. Pressure balance would be achieved if the electron 
temperature were 507o higher on side 2 than on side 1, but this is 
unlikely. 

We conclude that the discontinuity which passed Pioneer 6 at 

0109 UT on January 19 1966, was a reverse slow shock. 

A Forward Slow Shock . Figure 3 and 4 shows «o discontinuity 

which Ivaa the signature of a forward slow shock, namely , an increase 

in n, T and V and a decrease in B. It was seen at Pionaar 6 

atiulA^l li.’Bx dsv January 20, 1966.. 0 , The observed parameteEa 

and changes are consistent with ♦'he Rankine-Kugoniot equations, as 

indicated by the ics agreement between the observed values and a 

solution to the R-H equations found by Chao's procedure. The "predicted" 

values of are in good agreement with the observed values (Table 2). 

The shock normal > and speeds are given in Table 3, which also shows 

V and V , the nosrmalpflow seeadEv relative to able shOoku Note that 
* "2 

Vn t<i V <-ast , indicating that the fluid was decelerated in the shock 
”2 ’'l 

front. The Mach numbers, and computed from the numbers in 

Table 1, are sho%m in Table 4. They show the characteristics required 
for a slow shock. Computed estimates of the electron temperatures are 
also given in Table 4; it can only be said that they have the correct 
order of magnitude. 

All the available evidence suggests that the discontinuity In 
Figures 3 ts a ferwnrdr^Iow afaook 
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V. Electron Temperature 

Many of the discontinuities which we found in our search for 
shocks did not satisfy some of the necessary conditions for shocks, and 
are thus probably tangential dlacontlnuitlat or rotational discontinuities. 
The pressure is constant across both tangential and rotational 
discontinuities. This condition can be used to solve for the electron 
temperature. 

Measurements of the electron temperature near 1 AU, by Montgomery 
et al. (1968), showed that remains in the rarge of 10^ to 1.5 x lO^^K 
even though the proton temperature varies from 3 x 10^'^K to 10^°K. 

These "direct" measurments are based on a fit to a segment of the 
Maxwelliaji distribution. At low energies the spectrum is contaminated 
by photoelectrons from the spacecraft and there is a non -Maxwellian 
tail at high energies, so there is the possibility of some uncertainty 
in the direct measurements. 

An indirect method of measuring electron temperature, beised on the 
pressure batlance at tangential discontinuities, was introduced by 
Burlaga (I968), who found 10 ^°K. Although this is consistent with 

the direct measurements it was based on preliminary proton temperatures 
which were near the threshhold of the detector. Ogilvie and Ness (1969) 

introduced another indirect method, which gives however, they 

did not give an electron temperature. 

Since it is desirable to confirm the direct measurements by a method 
which is not affected by the photo sheath or non-Maxwellian tail, we have 
extended the method of Burlaga (1968), using the discontinuities that we 
have found in the search for shocks. 
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Pressure bslence gives the following condition 

We do not hsve roessurements of alpha particles , but for an estimate of 

the average electron temperature it is sufficient to use the average 

values, n /n - .045 and T /T ■2.5 (Hundhausen, 1968), since the 

contribution of the alphas to the pressure is relatively small. The 

measured proton tempera*'ure is that along the earth-sun line, T^, 

%rhile the appropriate temperature in (1), T^, generally differs from 

this value because of the proton anisotropy, T^^/T^«1.9 (Hundhausen, 

1968). Since the discontinuity surfaces are usually along the spiral 

direction (^urlaga, 1968), we may set T^«1.5Tj^« This leaves 2 unknoims, 

T , and T As a zeroth approximation we assume that T ,aT _. This 

is consistent with the observations of Montgomery et al. (1968) which 

show that the shape of the Maxwellian part of the spectrum does not 

change appreciably with time. With T *T , (1) then gives 

el e 

T - •0'44 <"2”g ' (2) 

® - t>g .OSTSCn^-n^) 

where w is the most probable thermal speed wV^2kT/m in km/sec, n is in 
cm , and B is in units of lO gauss. Clearly, (2) is applicable 
ot\i|f •nfiytWet'iretyli «a4nlngfu4scalyilfinaheiugfei:tainty 

ift^T^ils^^felWlvel^ilfcall'3 T^ is relativel) small. 

The discontinuities for the madaulandoti of T were selected from 

e 

the total set described in Section II by eliminating those which were 
identified as shocks and those for which the uncertainty in was 
large as a result of the accumulation of errors in (2). Our crit ion 
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for a large uncertainty is that determined from (2) with |a b|/B*.1, 
jAn|/n-,2, and 1at|/T-. 3> is greater than .5 x 105oK. T\d s left 31 
dis'^ontinui ties from which we obtained the distribution of electron 


temperatures shown in Figure 4a. The most probable value and the 

average value are both 1.8 x lO^^K, The width of the distribution is 

probably statistical rather than real and due to the propeigation of 

errors in the man^ algebraic operations in (2), The half width at 
5o 

half majdmum 10 K gives a measure of the spread in T^. Since there 

4 

are 31 measurements of T , the uncertainty of the average is +2 x 10 ©k, 

e “ 

Thus using the formal criterion |at|<5x lO^^K we obtain T^ - 

(1.8 +.2) X lO^^K. A more restrictive criterion for AT gives a somewhat 

smaller value of T . 

e 

The above method for estimating T^ has 2 difficulties: 1) it 

assumes that T^^^ = very sensitive to errors in 

n, w, and B so that only a fraction of the discontinuities can be 
used to cailculate Tg. As an alternative method, we assumed that T^^ = 
lO^^K or 2 X lO^^K and solved (1) for T^^ each case. The resulting 
distributions of T^^ shown in Figure 4b and 4c, In both cases, the 
median- i-.' value of (calculated) « T^^^ (assumed), consistent 
with our earlier assumption that the electron temperature usueilly does 
not change across a discontinuity. As in Figure 4a» the spread in 
the distr^outions is probably due to the propagation of errors. In 
any case, the symmetry of the distributions in Figure 4b euid 4c implies 
that T^ does not preferentially increase or decrease across the 


12 


discontinuities and the width of the distributions gives an upper 
limit. 10 °K, for the change in the temperature across a 

discontinuity. 
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V . Sonnnary 

We have searched through Pioneer 6 data flor 5 ' • •!.. 

reverse fast shocks, and forward and reverse slow shocks. No reverse 
fast shocks were found. A reverse slow shock was Identified; such a 
shock has not been reported heretofore. A forward slow sh'Ock waa . 
also Identified and v-was found to be similar to the twcxjfomrwL early by 
Chao and Olbert (1970). 

Numerous (50) other discontinuities were found, most of which are 
probably tangential since all of the plasma parameters (n, V, T) changed 
across most of them. In partlci lar, changes In B were essentially 
always accompanied by changes In n, which Indicates that discontinuities 
of the type suggested by Ivanov (with n^/n^ * 1) seldom If 

ever occur In the solar wind. More precisely. In a sample of discontinuities 
In B from 4 months of data from Pioneer 6, <57o were associated with 
no density change, l.e., <57o were the type of dlscontlnulttes predicted 
by Ivanov. 

Assuming that the non-shock discontinuities were either tangential 
or rotational and assuming no change in the electron temperature 
across the discontinuities, one can use the pressure balance condition 
to calculate T^. Using those discontinuities for which we could compute 
reasonably accurate (+507o> values of we obtained an average 

* (1.8+.2) xlO^^ from 31 of the discontinuities. This is in good 
agreement with the directly measurements values. 
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Figure 1 


Figure 2 


Figure 3 
Figure 4 


FIGURE CAPTIONS 

Distribution of discontinuities 

selected from Pioneer 6. The dip in the distribution 
of is due to the selection procedure which eliminates 

discontinuities with small changes in B. The n^^/n^ distri- 
bution indicates that disconclnuoufi changes in B are 
essentially always accompanied by changes in n. The 
lowest histogram shows that the temperature (thermal 
speed, w) tends to change also; it is not known to what 
extent the width of this distribution is determined by 
measurement uncertainties. 

Reverse slow shock of 19 January, 1966. Plasma and 
magnetic field measurements are shown together with best 
fits (solid lines) to the Rankine-Hugoniot equations, 
determined by Chao's procedure. 

Forward slow shock (See Figure 2). 

Computed electron temperatures, a) Distribution of 
temperatures computed from the pressure balance condition 
on the assumption that T^^ = and the requirement 

that AT^ < 5x10 °K. b) Distribution of temperatures 
T^ computed assuming that T^^ = 10^°K. c) Distribution 
of temperatures computed assuming that ** 2 x105ok, 


TABLE 1 


★ 

BEST-FIT VALUES 



Time 



Bit 


CVJ 

PQ 

^2T 

®2N 

®1 



V 

IT 


V 

2R 

V 

2T 

V 

2N 

’^1 


"1 

^2 

19 

Jan. 

0109 

4.4 

-2.0 

-0.2 

4.0 

-1.0 

1.0 

5.2 

4.2 

321 

15 

15 

310 

25 

22 

323 

311 

7.0 

9.0 

20 

Jan. 

1423 

-5.0 

6.0 

-4.0 

-4.5 

1.8 

- 3.6 

10.0 

6.0 

388 

13 

-10 

405 

-30 

4 

388 

4 o 6 

7.5 

10. 0 


0. i , 

0819 

3,^ 

9 . 7 

5 9 

'-*08 

<'.3 

•0 

11 .8 


^ . ■ 

\ 

3 

< .8 


8: 

U8 

36 3 

^ t 

- 


*Here B is in gauss, V is in km/ sec, and n is in cm . 


TABLE 2 


V^, Predicted and Observed 


Date 

''a 

^2R 

V 

V2t 

V 

2N 

km/sec 

January I9 

312 

311 

25 

23 

Theory 

0109 UT 

31U + 5 

310 + 5 

18 + 10 

5 ± 20 

Observed 

January 20 

4o6 

4o5 

-28 

4 

Theory 

1423 UT 

415 + 10 

4l4 +10 

-10 + 15 

0 + 20 

Observed 

October 27 

363 

348 

56 

87 

Theory 

0819 

359 + 10 

348 + 10 

55 + 10 

69 + 20 

Observed 



< c <c 


TABLE 3 


BASIC SHOCK PARAMETERS 


January I9, i960 January 20, I966 


V (km/sec) 

(212,1,135) 

(21*1,95,193) 

Vg (km/sec) 

251 

323 

n (RTN) 

(.85, .01, .5h) 

(.75, .30, .60) 


= 33°, jZis - 0° 

e = 37°, 0 = 22° 
s 

Vni*(Wsec) 

-29 

37 


-23 

27 

%,n 

42° 

o^ 

0 

0 


\ 


TABLE 4 

ADDITIONAL SHOCK PARAMETERS 


January 19, I966 


M 

A,1 

.9 

^A,2 

.8 

^lOWjl 

1.2 

^slow,2 

0.8 


«.4 

Te2(10^°K) 

».6 


Janusiry 20, I966 

.9 

.8 

1.3 

0.8 


« 11 


NUMBER 
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Figure 4 





